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We discuss friction-dominated vortex-string evolution using a new analytic model recently devel- 
oped by the authors. By treating the average string velocity, as well as the characteristic lengthscale, 
as dynamical variables, we can provide a quantitative picture of the complete evolution of a vortex- 
string network. Previously known scaling laws are confirmed, and new quantitative predictions 
regarding loop production and evolution are made. 
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I. INTRODUCTION 



The concept of symmetry breaking plays a crucial role in modern physics, and one of its most interesting conse- 
quences is the formation of topological defects. These defects have beeix-observed and studiedixi a wide variety of co 
densed matter contexts, including metal crystallization^, liquid crystalsoLl, superfluid heliumalj and superconductor, 
^ In models where they are allowed, defects will form whenever the rate of the phase transition is fast relative to the 
J> . scale of the system size (in other words, a 'quench'). 

CO ' On the other hand, they arg also believed to have formed in the early universe, and they can play an extrernely 
important part in its evolutionQ. In this context, the conditions for their forariation were first established by KibbleQ — 
' except for some subtleties in the case of the breaking of a gauge symmetryO, they are entirely analogous. 



The scaling evolution c^VjOjtex-string networks has been extensively studied analytically in both condensed matter 
and cosmological settingsE2rll3, but using rather different methods. This difference is perhaps understandable given the 
I extremity of these two physical regimes, but it may not be necessary. Condensed matter descriptions tend to focus on a 
"Y^i ' coarse-grained order parameter (/>, providing a low-level picture of defect motion by estimating energy dissipation rates. 

On the other hand, high energy physicists take an 'idealized' one-dimensional view of string dynamics by integrating 
^ , out the radial degrees of freedom (in the Higgs (j) and other fields) to obtain a low-energy effective action — the Nambu 
I ' action. The resulting relativistic equations of motion can then be averaged to describe the large-scale evolution of the 
' O \ string network. Naturally one should also account for energy loss mechanisms, such as loop production — something 
^ ' that is not done in condensed-matter contexts. _ 
Q \ The analytical study of casmological string netwocks was started one decade ago by Kibble's 'one-scale' modelEHl 
• • . (later modified by BennettB and Albrecht & TurokEJ). There it was assumed that the evolution of the long-string 
. ^ ' network could be described by a single lengthscale, which is usually called the 'correlation length'. One then supposes 
^ , that a scaling solution exists at late times and ends up showing that such a solution will in fact exist and be stable 
■ subject to conditions on the loop production mechanisms. One of the caveats of this model is that it is only valid in 
the latter evolutionary stages when strings arc moving rclativistically. However, at early times the string dynamics 
is friction-dominated, due to particle-string scattering. Of course, this is also the relevant case in condensed matter 
physics. |--.|--| 
A model of string evolution including the effects of frictional forces has been recently proposed by the authorsE^llj. 
It is a simple generalization of the 'one-scale' model in which the average string rms velocity becomes a dynamical 
variable. This simple-madel provides the first quantitative description of the complete evolution of a string network 
in the early universecjllj. In this paper, we will use it to study string evolution at constant temperature, which is 
relevant in condensed matter contexts. The main advantage of this approach is that one can easily account for the 
production and evolution of string loops. Although it should be seen as the basis for further work, the model is 
already predictive enough to be testable in laboratory experiments. 
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The structure of this paper is as follows. In the next section, after a review of string dynamics, the evolution 
equations for the 'characteristic lengthscale' and the average velocity of the long string network and each individual 
loop are derived and justified. The cases of strings arising from the breaking of gauge and global symmetries are 
both considered. The validity of these 'averaged' evolution equations is then tested against a simple — but physically 
relevant — loop solution in section 3. In section 4 we describe string evolution at constant temperature. The case 
where the dynamics is friction-dominated is directly applicable to a whole range of situations in condensed-matter 
systems. It will be shown that this model can easily reproduce some well-known results, notably the L cx t^/"^ law. 
Furthermore, new quantitative predictions regarding loop evolution are made. We also discuss the opposite limit, 
where strings are initially 'free' and find that in that case a network quickly evolves towards the previous regime. We 
note that the evolution equations for the string network and the string loops are relatively straightforward, but the 
intervening derivation is not. First-time readers may wish to skip straight to the results in section IV. 



II. THE AVERAGED EVOLUTION EQUATIONS 



A. Basics of string dynamics 

The usual condensed matter approach to vortex dynamics is based on a 'coarse-grained' complex scalar field In 
quantum field theory, for example, one can consider the abelian-Higgs model, which is a relativistic generalization of 
the Ginzburg-Landau theory of superconductivity. It is also of interest to consider the global version of this, that is 
the Goldstone model. I-I 

In high-energy physics it also proves to be convenient to adopt a one-dimensional view of string dynamics (seeQ 
for a detailed discussion of the issues below). In this description a string sweeps out a two-dimensional surface (the 
worldsheet) which can be described by two 'worldsheet coordinates' — one is time-like and can in fact be identified 
with the background time (which we will denote by i), while the other is space-like and simply labels points along 
the string (we will call it a). 

This one-dimensional description is essentially acheived by integrating over the radial modes of the vortex solution 
on the assumption that the scale of perturbations along the string is much larger than its width — thereby obtaining 
a low-energy effective action. For the case of a gauge (global) string, one thereby obtains the Nambu (Kalb-Ramond) 
action from the abelian-Higgs (Goldstone) model. Note that in the global case there are long range forces between 
the stringfi. in general there can be both external and self-field contributions to this. However, it can be shown 
(again, seeQ for a detailed discussion) that the later can be 'rcnormalized out', yielding the well-known logarithmically 
divergent string energy per unit length /i (which is a constant for gauge strings). By varying these actions it is then 
straightforward to obtain the string equations of motion. 

There is, however, a crucial ingredient for string evolution missiujg^ Since strings move through a background fluid, 
their motion is retarded by particle scattering. Vilenkin has shownEZi that this effect can be described by a frictional 
force per unit length that can be written 

Ft = -^7v, (2.1) 

where v is the string velocity, 7 is the Lorentz factor and F is a constant damping coefficient, that can be written as 
the square of a characteristic propagation speed (which need not necessarily be the speed of light) times a 'friction 
timescale' rt, whose explicit value depends on the type of symmetry involved. For a gauge string, the raain contribution 
comes from Aharonov-Bohm scatteringl3, while in the global case it comes from Everett scatteringij. For example, 
if the background fluid is a perfect gas, we have 

( 2-Kh (ksT^f 



I ¥^te^l-(f)ln^(T^'5) Global 



where T is the background temperature and /3 is a numerical factor related to the number of particle species interacting 
with the string (strictly speaking, its value is slightly different in the two cases, but a common symbol will be used 
for simplicity) . It should also be noted that the Everett scattering formula is only valid when the particle wavelength 
is much larger than the string thickness 6. Alternatively, F can be written as the characteristic propagation speed 
times a 'friction lengthscale' £f. 



For stritLg motion in a flat background, the string equations of motion with the frictional force (2.1) can then be 
written asllj 



2 



(2.3) 



0, 



(2.4) 



where the dimensionless parameter e (which can be interpreted as a 'coordinate energy per unit length') is defined by 

e'^ , (2.5) 

and dots and primes respectively denote time and space derivatives. This form of the evolution equations proves to be 
particularly useful because dissipation is naturally incorporated in the decay of the coordinate energy density e, while 
preserving the gauge conditions. Note that while this is a truly relativistic formalism, it is straightforward to obtain 
the non-relativistic l imit tha t will be adequate to condensed matter contexts where the dynamics is friction-dominated. 
In this case, Eqns. (2.3-^) reduce to 



X 

f 



-[x'A(x'Ax")] , 



(2.6) 



and one can recognize the right-hand side as the friction force term (which is dominant in this limit), e.g. on a 
superfluid vortexEJ j_. 

Incidentally, it has been showrciJ that a global string will behave as a superfluid vortex if it is introduced in a 
homogeneous background (with a density ph , say) — physically, this corresponds to giving it angular momentum. The 
interaction between this background and the string originates an additional force, known as the (relativistic) Magnus 
force, and (2.3) becomes 



1 - 



(2.7) 



1/2 

where T' oc ; note that (2.4) remains unchanged 



B. Lengthscale evolution 



We can now proceed to average the string equations of motion to describe the large-scale evolution of the string 
network. We therefore define the total string energy and the average rms string velocity to be 



E ^ p 



sda . 



(2., 



2 _ /^.^2\ _ J^'^eda 



V = X = 



Jeda 



(2.9) 



Differentiating (2.8) and using (2.4) and (2.9), we see that the total string energy density will obey the following 



evolution equation: 



dp 
'dt 



0. 



(2.10) 



This includes both long strings and loops which have a low probability of interacting with other strings before 
decaying. We shall study the evolution of the long-string network on the assumption that it can be characterized by 
a single lengthscale L; this can be interpreted as the inter-string distance or the 'correlation length'. Strings larger 
than L will be called long or 'infinite'; otherwise they will be called loops. For Brownian long strings, we can define 
the 'correlation length' L in terms of the network densitjf] p^o as 



^Throughout this paper the subscript 'oo' refers to properties of the long ('infinite') string network. 
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Poo^^- (2.11) 

Following KibbleQ, the rate of loop production from long-string collisions can be estimated as follows. Conceptually, 
we divide the network into a collection of segments of length L, each in a volume . Consider another segment of 
length I moving with a velocity Voo] the probability of it encountering one of the other segments within a time 5t is 
approximately IvooSt/L'^. Consistently with our 'one-scale' assumption, we then assume that the probability of such 
an intersection creating a loop of length in the range I to I + dl will be given by a scale- invariant function w (l/L). 
The rate of energy loss into loops is then given by 

dpoO \ Voc f f i\ £ d£ ^ poo /O 1 r,^ 

= Poo^ / W - -— = CToo^. (2.12) 



dt /to loops L J \LJ LL - L 



where the loop 'choppi ng' e fficiency c is assumed to be constant. By subtracting the loop energy losses ( ^.12| ) from 
(2.10) and then using ( ^.11 ), we obtain the overall evolution equation for the characteristic lengthscale L, 



dT 7)^ 

2—^^L+~cVoo. (2.13) 
dt 1 



C. Loop evolution 



The main advantage of the present approach is that we can also study the evolution of the loop density and 
distribution, which have so far been neglected in the condensed matter literature. The traditional approach in 
cosmology is to define ni{£,t)d£ to be the number density of loops with length in the range {£,£ + d£) at time t\ the 
corresponding loop energy density distribution is 

pei£,t)d£ = fi£ni{£,t)d£. (2.14) 

Note that the total loop energy density is 

Po= [ Pii£,t)d£, (2.15) 



and p — Poo + Po0 From our assumptions on the loop production rate ( 2.12| ) it is then easy to see that 



dpe , vj Voo£ ( £ \ 



dt ' V"' \L 

However, note that this equation is 'static', in the sense that it does not include loop decay mechanisms. 

Instead, we start by using our analytic model to describe the evolution of each individual loop. Knowing the energy 
density transferred from long strings into loops and estimating their sizes at formation (see below) , one can determine 
the energy density in loops and other relevant quantities at all times. This formalism does not allow for possible loop 
reconnections or self-intersections but these should not be important in a friction-dominated regime. 

The physical size of a loop is simply given by 

£ = I eda, (2.17) 

J loop 

and its time evolution is easily found to be 



■^Throughout this paper the subscript 'o' refers to properties of the entire loop population, while refers to the loops with 
length in the range (i, I -\- dl). 
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Now, we will assume that loop production is 'monochromatic', ie that loops formed at a time tp have an initial 
length 



e{tp) = a{tp) L{tp) , (2.19) 



where a is a parameter of order unity. 



With this ansatz the scale-invariant loop production function w (see (2.12)) becomes 



yjli)=^s(^-a] , (2.20) 



L J a \L 



and the rate of energy loss into loops becomes 



dpoo \ ^ ^g^oo 

'^^ ) to loops 



Hence the energy density converted into loops from time t to t + is 



dpo{t) = ^ic^dt; (2.22) 
this corresponds to a fraction 

7-T = c— dt, (2.23) 

Poo(i) L 

of the energy density in the form of long strings at time t. Then using our ansatz ( ^.20 ), the corresponding number 
of loops produced in a volume V is 

dN{t) = -'^Vdt- (2.24) 
hence the ratio of the energy densities in loops and long strings at time t is (neglecting the initial loop population at 

t = Q 

where tc is the moment of the network formation and £(i, t') is the length at time t of loops produced at time t' . 

We can therefore numerically (and, in some simple limit cases, analytically) determine the loop density at all times. 
This generalized 'one-scale' model can therefore provide a complete description of a string network. 



D. Velocity evolution 

We now consider the evolution of the average string velocity v. A non-rclativistic equation can be easily obtained: 
it is just Newton's law, 



P dv n u 



This merely states that curvature accelerates the strings while friction slows them down. On dimensional grounds, 
the force per unit length due to curvature should be /i over the curvature radius R. The form of the damping force 
can be found similarly. 

A relativistic generalization of the velocity evolution equation ( 2.26 ) can be obtained more rigorously by differen- 
tiating dJ): 

1 dv / / k v\ 

This is exact up to second-order terms. In the curvature term, we have introduced R via the definition of the 
curvature radius vector, 
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where u is a unit vector and s is the physical length along the string (related to the coordinate length a by ds 

1 /2 

jx'jdCT = (l ~ x^/c^) eda). The dimensionless parameter k is defined by 



~2 



((l--)(i.u))^fc« 1--) . (2.29) 



Note that in the case of long strings, our 'one-scale' assumption implies that the curvature radius coincides with 
the correlation length, R = L;on the other hand, for a loop of size £ we should have £ « 2wR. 

The parameter k is 'phenomenological', and is related to the presence of small-scale 'wiggles' on strings: on a 
perfectly smooth string, u a*«i x will be parallel so fc = 1 (up to a second-order term as above); however this need 
not be so for a wiggly stringla. On the other hand, one can show that in flat spacetime (with no friction) fc = 0. 

In the case of long strings, we should expect fc = 1 since the correlation length L is much larger than the 'friction 
length' £f. Now, consider a particular string loop. While it is large compared with the friction lenghtscale, taking 
fc « 1 should again be a good assumption. In the opposite limit friction is no longer effective in damping its motion, 
and it moves as if it was in flat spacetime. Thus one requires that fc — > as i? — > 0. In particular, demanding that 
their limiting velocity be vj/c^ — 1/2 leads to the requirement that k (x R as R ^ 0. With these requirements in 
mind, and after comparing with the 'microscopic' (ie, unaveraged) evolution of some simple solutions (to be described 
in the next section) one arrives at the following ansatz: 

k=l]lR til (2.30) 

1^ V2 ' il ^ ^ 

where x is a numerical coefficientii of order one. Recall that the physical loop length is approximately £ = 2ttR, while 
for the long-string network L = R. 

Equations (2.13), (2. IS) and ( p. 27 ) form the basis of our generalized 'one-scale' model, which we will now proceed 
to apply. 



III. 'AVERAGED' VERSUS 'MICROSCOPIC EVOLUTION 



In order to check the validity of our 'averaged' evolution model, and in particular our ansatz for fc, we will test it 
against a simple loop solution. 

Consider a circular loop in a flat background with a constant friction timescale — that is, a condensed-niatter-like 
situation. We can describe the loop trajectory simply by 



X = r(i)(sin6',cos6l,0) , 9 G [0,2tt] 



Then equations (p.3|,2.4) reduce to 



f2 



1 



Note that the physical ('invariant') loop radius is i? = r j \J\ ~ j (p- , obeying 

R 



-4^ 



also the 'microscopic' velocity is w = — f and obeys 

4i' = ( 1 - 



1 V 
r ^ f 



On the other hand, our averaged evolution equations ( 2.18 , 2.27 ) take the form 



dR -TP- 

-di=-^T 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



(3.5) 
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1 dv _ f /k{R) V 



(3.6) 



Notice the similarity between the two approaches. Loops with size much larger than a 'friction lengthscale' that can 
be defined in the obvious way from the friction timescale will be overdamped, with the velocity being approximately 
given by 

v^^. (3.7) 

In this case the two sets of evolution equations actually coincide — hence justifying our k = 1 ansatz for large R. 
As the loop gains velocity r and R become significantly different and this equivalence ceases to be valid. When R 
becomes much smaller than £{, the loop still looses energy due to friction, but this is no longer effective in damping its 
motion — the loop now begins to oscillate relativistically. In particular, over one 'period' v oscillates between and 1. 
But we know that the averaged velocity should (in the limit) be iP' jc? = 1/2; this is the physical reason why we need 
a 'phenomenological' varying k on s mall scales. As we mentioned previously, this requirement fixes the behaviour of 



k on small scales to be as shown in (2.30). The remaining question is then how to match the two regimes. 

First of all, we need a clear idea of when (and where) the transition occurs. A good guess would be the moment 
of the 'first collapse', ie, the moment when we firstpiave v — c. In fact, this turns out to be a well-defined event. As 
was first pointed out by Garriga and SakellariadouEj (and can be easily seen by analytical or numerical study of the 



equation of motion (3.2)), circular loops with initial radius much larger than the friction length always reach v = c 



for the first time when 



0.5691. (3.8) 



Note that ^ if is the physically relevant case for string dynamics in condensed matter contexts (recall that the 
dynamics in that case is always friction-dominated). Also note that because of friction, all loops will rapidly become 
(almost) circular. 

After numerically comparing the averaged and microscopic evolution equations, we find that the simplest possibility, 

k = l^ln til (3-9) 

(see figure |^) provides the best answer (see figures In particular, this turns out to be significantly better than 

assuming smoother (and slower) transitions between the two regimes. As can be readily seen, this ansatz provides a 
very good fit, considering the lack of parameters available. _ 

In passing, it is worth pointing out that one can also easily calculate the loop lifetimec3. In the relativistic regime, 
the R evolution equation can be written 

f = --' (3-10) 
at T[ 

so we can immediately estimate that the loop will disappear in a time tdec ~ 2Tf after its first collapse. 

Having thus established, in a simple but physically relevant case, the validity of our generalized 'one-scale model', 
and in particular of the ansatz for fc, we now proceed to apply it to the study of string evolution in condensed matter 
contexts. 

IV. STRING EVOLUTION AT CONSTANT TEMPERATURE 
A. The condensed-matter context 

1. Network scaling 

As we already noted, in this case the dynamics is always dominated by friction. This means that the 'correlation 
length' L should always be larger than the (constant) friction length, so we can take fc = 1. Then the evolution 
equations can be approximated by 
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the friction lengthscale and the string energy per unit length being respectively 



. 1 s Gauge , , 

^'-\ slnff) Global (^-^^ 



and 



r r2 Gauge 
'^=|rhn(f) Global (4.4) 

where is the temperature at which the strings form and s is a constant. We then find the following late-time 
asymptotic behaviour: 

^/TT^(^^)^/^ (4.5) 

v = ^. (4.6) 
Note that in both cases the asymptotic behaviour of the long-string density is 



{l + c)hscH ' 



(4.7) 



the extra logarithmic dependencies in the global case cancel out. It should be. emphasized that in condensed-matter 
analyses one does not consider loop formation, although there is experimentalEI and computational evidence for them. 
Our results show that loop formation can play an important evolutionary role. The asymptotic ratio of the loop 
production and friction terms is a constant, which is precisely c — which in this way acquires a clearer physical 
meaning. As expected, increasing c (or including loop losses in the first place) leads to a lower network scaling density 
and a smaller average velocity v] furthermore, the approach to the scaling regime is also faster. 

Figure ^ illustrates the relaxation to scaling of a network of gauge strings for a particular set of initial conditions, 
with and without loop production. The differences between the two cases are clearly visible. Notice that in the 
gauge case the temperature only enters the scaling solution in the prefactor s; hence the plots of figure ^ (where the 
dimensionless variable L/ s is used) are 'universal curves', that is, hold for all relevant temperatures. 

In the global case, the logarithmic dependence of the friction lengthscale gives rise to an additional logarithmic 
dependence of the scaling solution on temperature. One can define d via L/5 = L/ sd] in the Goldstone model we 
then have 

The approach to scaling of a network of global strings is shown in figure ^ for two different temperatures. Note the 
enhancement of loop production at the early stages, since the string velocity is high; correspondingly, there is a fast 
growth of the correlation length. Comparing with the gauge case (see figure ^ , one finds that the effect of the extra 
logarithmic terms is significant for at least three orders of magnitude in time. 



2. Loop populations 



Now let us consider the loop populations. In the traditional (that is, 'static') approach one would simply calculate 
the loop distribution pe{£,t) in the scaling regime — which we assume starts at some time ts- The corresponding 
evolution equation is then 



dpi v'j , ^ . , ^ 

' " 5MT5^ T ■ (4-9) 



dt ' T^' ""^L^^^L 

Since the loop size and velocity are initially (that is, while the loops are overcritically damped) slowly-varying quanti- 
ties, this approach will only be accurate at the early stages of evolution of each lengthscale £. Afterwards £ and v will 
vary quickly (the averaged v will be a constant later), and in this situation this approach is no longer useful. Note 
that since friction is dominating the dynamics, irregularities are quickly erased and all loops become circular — hence 
the discussion of section 3 is particularly relevant here. 

Consistent with what we found there, we will use the following simplifying ansatz 

(where £ — 2ttR); naturally in this approach only the first case is relevant. Furthermore, the loop production parameter 
a (defined is ( 2.19[ )) should be a constant (we are considering a scaling regime) of order unity. 



After some algebra, we find the following solution 

p,^e, t) = I H^'^^^t iw - Jl^)] aLs<£< aL ^^^^^^ 

' l_ otherwise 

where we have neglected any loops present at time tg. Note that in the global string case these expressions are only 
approximate, since there are additional logarithmic dependencies in the energy per unit length fi and the friction 
length £f. 

However, we can do more than that. To begin with, we can determine the loop lifetime. During most of it the 



loop will have a length greater than the friction lenghtscale, so we can approximate their velocity by (4.10). Then its 
length will vary according to 



dt 

so its lifetime is 



d£ , , 



t{£) _ W ^ 

Tf ~ 87r2 Uf 



(4.13) 

Of course a more accurate result can be obtained numerically (see figure ^ ; for large enough loops we find 



2 



n ~ \£f 

showing our analytical estimate to be correct to within a fraction of one percent 



M^l ,^-0.0126, (4.14) 



On the other hand, we can simply use (2.25) to determine the ratio of the energy densities in the form of loops and 



long strings. In this case, taking a ~ 1, it simplifies to 

m = ^ = cL\t) f "-f^£{t, t')dt' . (4.15) 

In fact, an approximate analytic solution for g can be found in the scaling regime. We will assume that such a 
regime starts at a time ts and neglect any existing loops at ts- The contribution of each loop to the loop density 
will only be significant while it overdamped. In this regime we can assume that its length is approximately constant, 
£{t,t') ~ £{t') = L{t') (hence we will be deriving an overestimate of g). O n the other hand, the end of this regime 



approximately coincides with the moment of loop decay, which according to ( 4.14 ) for a loop formed at time t' happens 
at 

t = C^t' ::^[1 + L{l+C)]t' , (4.16) 

for a loop formed at time t' . Then we have 

^ ' ' 10 otherwise ^ ' 
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and there are two corresponding cases to consider for the integral (4.15) (see figure 0), according to whether or not 
loops have decayed at the time in question, 

git)^cL^{t) f "^df, (4.18) 



where 



We then obtain 



= <llt fi (4.19) 



g.caUng = { ^ ^) ^ ^ ^ (4.20) 

[id t> ^-H, 

This is easily confirmed numerically. Figure ^ displays the evolution of the ratio of the loop and long string densities 
corresponding to the gauge string network of figure ^ In this case c = 1, so our analytical estimate is Qscaiing = 0.0126, 
while numerically we find Qscaiing ^ 0.0089. Notice that with our initial conditions the ratio of the energy densities 
in loops and long strings is larger while the network is approaching scaling than in the scaling regime. This is simply 
due to the fact that loops with initial sizes not much larger than the friction lengthscale live longer than predicted by 



(4.14) (as can be seen in figure |6|). 

Naturally, since we can find out how each loop evolves, it is always possible to find out how many loops there are 
at any time on any given length interval. Note that in this context reconnections and self-intersections should be 
negligible, at least when the loops are large compared to If (because they are non-relativistic) — which is when their 
contribution to the loop density is significant. 

Furthermore, using the results from our discussion in section 2C, we easily find that the fraction of the energy 
density in long strings at time t that is converted into loops in the time interval At — Tf is asymptotically 

/n- (4.21) 

^ ' 1 + 5 t ' ^ ' 

note that this is a very small number. Furthermore, note that as we increase c there is a limiting fraction. Also, the 
number of loops formed in a time interval Ai = Tf and volume (.f is 

^ ; (4.22) 



(1 + c)-V2 V t 

in this case, Urt goes to zero as c grows — this is because L and hence the loop size at formation grow with c. 



3. Discussion 



The t^/^ scaling law for the characteristic lengthscale is a well-known result in the theory of phase ordering (that is, 
the growth of 'order' — as measured by some correlation length — by domain coarsening when a system is 'quenched' 
from a homogeneous phase into a brokeitsymmetry phase) with a non-conserved order parameter. In this context it 
is usually called the Lifshitz-Allen-Cahnc3 growth law, and it is widely supported by simulations and experiments 
(see also ref. |2|for a recent review). In the usual approach, one sets up a continuum description in terms of a coarse- 
grained order parameter and then assumes a scaling hypothesis^ that is, that at late times there is a single lengthscale 
such that the domain structure is time-independent (in a statistical sense) when all lengths are rescaled by it. The 
growthJaw is usually derived by studyijHs; the dynamics of the defects in (see, for example, Section 3 in Bray's 
reviewlla) . A recent alternative approached proceeds instead by comparing the global rate of energy change due to the 
energy dissipation to the local evolution of the order parameter; with the scaling hypothesis, the time-dependence of 
the lengthscale can be determined self-consistently. 

In particular, the L oc t^^^ law has been experimentally confirmed for the evolution of a string network in a 
nematic liquid crystal (roughly speaking, a liquid made of rod-like molecules) — eg, see ref. ^ where, as mentioned, 
loop formation and decay have been seen. 

The V oc L~^lnL scaling lawris also known in hydrodynamical contexts. Furthermore, it has been shown that it 
holds for superfluid vortex-ring£il in the context of a modified Goldstone model (in a way described in the previous 
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section). Hence the above result seems to indicate that a global string network at constant temperature asymptotically 
behaves as if it was made of loops of size L. 

As we hopefully made clear, the above discussion of the loop lifetime and density is entirely new in this context. It is 
hoped that these results, obtained in a fairly simple way in this approach will stimulate condensed matter theoreticians 
an experimentalists... 

Finally, in the c ase of supcrfluid vortices, despite the additional Magnus force term, the evolution equations also have 
the form (4.1,4.2). In this case, however, the physical meaning of the friction lengthscale is not clear. Furthermore, 
it is also not clear how one can describe the effect of the Magnus force on the evolution of the network. We hope to 
address these issues at a later stage. 

Therefore, a model aimed at describing the evolution of cosmic stringslijllj can be straightforwardly applied to rather 
more 'down-to-earth' situations, with some advantages over previously existing analytical approaches (in particular, 
allowing a precise treatment of loop production). This is rather remarkable, considering the difference in the energy 
scales involved. 



B. The 'relativistic' regime 

As a matter of completeness as well as mathematical curiosity, we now consider the evolution of a string network in 
flat space with a constant friction lengthscale when the initial conditions are such that the correlation length is much 
smaller than the friction lengthscale. This is therefore the opposite regime to that usually observed in condensed 
matter. 

The evolution equations will now be 

dL 

2-^-L + 5., (4.23) 



dv 
'dt 



fk{L) V 



(4.24) 



where we will assume that our ansatz for k for string loops also holds for long strings (see ref. 16 for a more complete 
discussion of this point), that is 



k = 



1' f >X 

J_ i ^ < Y 

f t ' t., ^ ^ 



(4.25) 



In the regime where i? <C the u-equation is independent of L, so its particularly easy to find the scaling regime 



L 



if 



c exp 



- 2V2c 



(4.26) 



Hence L grows exponentially fast and quickly 'catches up' with if, in other words, a network starting in the 'free' 
regime rapidly evolves to the usual friction-dominated regime. Figure |^ displays the evolution of a gauge string 
network in the free regime, and in particular the transition to the damped regime. Note that if loop production is 
allowed, this fast growth of the correlation length will obviously mean that an extremely large number of loops is 
produced. In this case, the energy density in loops actually exceeds the energy in long strings — in the above plot, by a 
factor of 10^. A word of caution is however needed here. In this case, loop reconnections onto the long string network 
should play an important role. However, since we still get an exponential growth if loop production is switched off 
(c = 0) — although the growth rate of L is obviously much larger for c ^ — our results should at least be qualitatively 
correct. 
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V. CONCLUSIONS 



In this paper we_.]iave discussed the apphcabihty in condensed matter contexts of a recently-developed model 
of string evolutionEjO, where a 'characteristic lengthscale' and the average velocity are the dynamical variables. 
This has allowed us to properly describe friction-dominated string dynamics, hence providing the first complete and 
fully quantitative study of string networks and their corresponding loop populations in condensed matter (as well 
as cosmologicaltj) contexts. The fact that these results can be obtained in a model initially aimed at describing 
cosmic string evolution is, of course, a manifestation of the universality of symmetry breaking and defect formation 
phenomena, but it also lends weight to the validity of this approach because these cosmological models have been 
extensively tested numerically. 

We have confirmed two previously known condensed matter scaling laws, while noting the significant effects of 
loop production (which is not considered in standard approaches). We also presented the first (analytic or otherwise) 
discussion of the evolution of the overall loop density, while also determining loop lifetimes. As we have already stated, 
the literature on string loops in phase ordering kinetics appears to be restricted to a few references to the observation 
of loop production and decay in laboratory experimentso. It is hoped that these more quantitative predictions will 
trigger further efforts from the condensed matter community. 
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The ansatz for k 
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FIG. 1. The ansatz for parameter k for string loops in a condensed-matter context as a function of the rescaled loop radius. 
The behaviour at large and small scales has a physical justification (see text); matching was done numerically. 
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(a) Loop radius (b) Loop velooity 
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FIG. 2. Comparing the 'microscopic' (solid lines) and 'averaged' (dashed) evolution equations for a circular loop in a con- 
densed-matter context. Length and time are in units of £i, and the time axis is with a logarithmic scale. Plot (a) depicts the 
log of the (rescaled) radius, while (b) depicts the loop velocity (in units of c). The lower graphs are close-ups of the upper ones. 
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(a) Loop radius (b) Loop velocity 
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(c) Loop radius (d) Loop velocity 




FIG. 3. Comparing the 'microscopic' (solid lines) and 'averaged' (dashed) evolution equations for a circular loop in a con- 
densed-matter context. Length and time are in units of it, and the time axis is with a linear scale. Plot (a) depicts the log of 
the (rescaled) radius, while (b) depicts the loop velocity (in units of c). The lower graphs are close-ups of the upper ones. 
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(a) (b) 




log (t/tc) log (t/tc) 

FIG. 4. The approach to scahng of a network of gauge strings at constant temperature for c = (sohd hnes) and 6=1 
(dotted hnes); k = 1 and the initial conditions are O.lLi = c±i = s, Vi — 0.1 c. The plots correspond to the evolution of L/s (a) 

and V (in units of c) (b) and the exponents of the power-law dependence of L (c) and v (d). The horizontal axis is labeled in 
orders of magnitude in time form the moment of string formation; in (a) an (b), the 'y' axis is also in a logarithmic scale. 



16 



(a) (b) 
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FIG. 5. The approach to scahng of a network of global strings at constant temperature for T/Tc = 0.1 (solid lines) and 
T/Tc = 0.01 (dotted lines); c = 1, = 1 and the initial conditions arc O.lLi = cU = s, Vi = 0.1c. The plots correspond to 
the evolution of the exponents of the power-law dependence of L (a) and v ih), v itself (in units of c) (c) and the ratio of the 
loop formation and friction terms (d) . The horizontal axis is labeled in orders of magnitude in time from the moment of string 
formation; in (c) an (d), the 'y' axis is also in a logarithmic scale. 



17 



Loop decay 
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FIG. 6. The evolution of loops of different sizes relative to the friction lengthscale. For large enough loops, the lifetime is 
proportional to the square of the length. 
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FIG. 7. The interval in t' giving a non-zero loop length contribution to the integral (4.15) for different times t 
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FIG. 8. The ratio of the loop and long string densities for a network of gauge strings with c - 
the corresponding case of figure ^ 



1. Initial conditions are as in 
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FIG. 9. The evolution of an initially free network at constant temperature, for c = (solid lines) and c = 1 (dotted lines). 
Plots represent the lengthscale relative to the friction length (a), velocity (in units of c) (b) and ratio of long-string and 

loop energy densities (c); all are log- log plots (time is in orders of magnitude after the network formation). Notice that the 
exponential growth in L (aided by the relativistic velocity) leads to a period of loop domination. 
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